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Abstract. We review recent developments in the theory of the Lya forest and their implications 
for the role of the forest as a test of cosmological models. Simulations predict a relatively tight 
correlation between the local Lya optical depth and the local gas or dark matter density. Statisti- 
cal properties of the transmitted flux can constrain the amplitude and shape of the matter power 
spectrum at high redshift, test the assumption of Gaussian initial conditions, and probe the evolu- 
tion of dark energy by measuring the Hubble parameter H(z). Simulations predict increased Lya 
absorption in the vicinity of galaxies, but observations show a Lya deficit within A, ~ 0,5h Mpc 
(comoving). We investigate idealized models of "winds" and find that they must eliminate neutral 
hydrogen out to comoving radii ~ 1 .5h Mpc to marginally explain the data. Winds of this magni- 
tude suppress the flux power spectrum by ~ 0. 1 dex but have little effect on the distribution function 
or threshold crossing frequency. In light of the stringent demands on winds, we consider the alter- 
native possibility that extended Lya emission from target galaxies replaces absorbed flux, but we 
conclude that this explanation is unlikely. Taking full advantage of the data coming from large tele- 
scopes and from the Sloan Digital Sky Survey will require more complete understanding of the 
galaxy proximity effect, careful attention to continuum determination, and more accurate numerical 
predictions, with the goal of reaching 5 — 10% precision on key cosmological quantities. 

PHYSICS OF THE FOREST 

The 1990s saw four epochal advances in our understanding of the Lya forest. Spectra 
of quasar pairs showed coherence over scales of a hundred kpc and more, implying 
large sizes and thus low densities for the absorbing structures yll2ll2l0]. Keck HIRES 
spectra of unprecedented resolution and signal-to-noise demonstrated the ubiquity of 
weakly fluctuating Lya absorption in the high redshift universe [5], and they revealed 
thepresence of metal lines associated with low column density hydrogen absorbers 
10,0]. Finally, and most directly relevant to this review, a combination of numerical 
simulations and related analytic models led to a compelling new physical picture of Lya 
forest absorption fiEE GH EEO HQ H- 

The basic numerical result is simple to summarize: given a cosmological scenario 
motivated by independent observations, 3-d simulations that incorporate gravity, gas 
dynamics, and photoionization by the UV background produce something very much 
like the observed Lya forest, an outcome that requires no ad hoc adjustments to the 
model. The top three rows of Figure [T] illustrate this point, showing, respectively, the 
observed Lya forest of the z = 3.62 quasar Q1422+23 1 , expanded views of four selected 
regions of this spectrum, and simulated spectra of the same length along four randomly 
selected lines of sight through a cosmological simulation. The simulation uses smoothed 



particle hydrodynamics (SPH) with 128 3 dark matter particles and 128 3 gas particles in 
aperiodic cube of comoving size 11.H1/j~ 1 Mpc (1422 km s _1 at z = 3). It assumes a 
ACDM model (inflationary cold dark matter with a cosmological constant), with Q, m = 
0.4, Q A = 0.6, h = 0.65, Cl b = 0.02/T 2 = 0.0473, inflationary spectral tilt n = 0.95, 
and a power spectrum normalization that corresponds to o% = 0.8 at z = 0. Spectra are 
extracted from the z = 3 simulation output using the methods of II 1211 . 

Studies of the Lya forest have traditionally focused on "lines" identified by a de- 
composition procedure and the "clouds" or "absorbers" that produce them. Hydrody- 
namic simulations show that typical marginally saturated absorption features at z ~ 3 
arise in filamentary structures, which are analogous to (but smaller in scale than) to- 
day's galaxy superclusters. However, in the simulations there is no sharp distinction 
between the "lines" and the "background," and one can also characterize the Lya for- 
est as "Gunn-Peterson" lfl6ll absorption produced by a smoothly fluctuating intergalactic 
medium dO, EE d. 

What makes the fluctuating IGM perspective especially powerful is the tight relation 
between the density and temperature of low density cosmic gas, T ~ Tq(p /p) a with 
a ~ 0.6, which emerges from the balance between photoionization heating and adiabatic 
cooling [20, 21]. The Lya optical depth of photoionized gas is % °c hri °= p 2 r _0 ' 7 r _1 , 
where r 7 accounts for the temperature dependence of the recombination rate and 
r is the rate at which neutral atoms are ionized by the cosmic UV background. The 
temperature-density relation then leads to the "Fluctuating Gunn-Peterson Approxima- 
tion" (FGPA), F = exp(— x) = exp[— A (p/p)P] , which relates the continuum normalized 
flux F to the local gas overdensity H17LI22LI18L 19]. The index /3 lies in the range 1.6 — 2 
depending on the thermal history of the gas, and the proportionality constant A is itself 
proportional to Qfyv'TQ Qn T~ l . Thus, one can think of a Lya forest spectrum as provid- 
ing a 1 -dimensional, non-linear map of the gas overdensity p/p along the line of sight. 
The map is smoothed by thermal broadening and distorted by peculiar velocities, but 
these effects are small enough to leave a tight correlation between Lya optical depth 
and local gas density even in redshift space fl7ll23ll . Furthermore, pressure gradients in 
the diffuse, photoionized IGM are usually weak compared to gravitational forces, so the 
gas overdensity traces the dark matter overdensity fairly well. 

The FGPA provides a valuable way of thinking about the information content of 
the Lya forest, and it can be a useful calculational tool if one has a way to produce 
realizations of cosmic density and velocity fields. The analytic models of the Lya 
forest by Ja |2jlli3l[l4j] essentially combine the FGPA with a log-normal or Zel'dovich 
approximation model for creating non-linear density and velocity fields. Alternatively, 
one can run an inexpensive, gravity-only N-body simulation and assume that the diffuse 
gas traces the dark matter (or add an approximate treatment of gas pressure 1 251 1261 l27lo . 
The dotted lines in the third row of Figure [l] show spectra extracted from a particle-mesh 
(PM) N-body simulation, run from the same initial conditions as the SPH simulation, 
illustrating both the effectiveness and the limitations of this approach. The two spectra 
trace each other over most of their length, with the largest breakdowns occurring in 
regions where shock heating has pushed the gas above the T p a temperature-density 
relation, which makes the optical depth in the SPH simulation lower than the FGPA 
predicts. 
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FIGURE 1. Simulating the Lya forest. The top panel shows a continuum normalized, Keck HIRES 
spectrum of the Lya forest region of the quasar Q1422+231 (z = 3.62), from [6J. The next row of panels 
shows blowups of four regions 1422 km s in length. The third row shows simulated spectra extracted 
along four random lines of sight at z = 3 from an SPH simulation (solid) and a PM simulation (dotted) of 
the ACDM model, with il„, = 0.4, h = 0.65, and ff 8 = 0.8 (at z = 0). The fourth and fifth rows show PM 
results for different cosmological parameter values: o% — 1.0 and D. m = 0.2, respectively. All simulated 
spectra are 1422 km s _1 in length. 



The N-body+FGPA technique offers a convenient way to investigate the response of 
the Lya forest to changes in cosmological parameters. Comparing the third and fourth 
rows of Figure [l] illustrates the effect of increasing the matter fluctuation amplitude by 
20%, to o%(z = 0) = 1.0, with the same Fourier phases in the initial conditions. The 
increased clustering of the high amplitude model is evident in the greater incidence of 
saturated absorption and in the merging of multiple small scale features into single larger 
scale features. The bottom row shows spectra from a model with Cg = 0.8 and the same 
linear power spectrum but a lower matter density, Q. m = 0.2. The reduction in Q, m lowers 
the value of the Hubble parameter H(z) at z = 3, and as a result features are more densely 
packed in redshift space, producing "choppier" spectra. In fact, our ll.Hl/f~ Mpc 
simulation cube is only 1024 km s _1 in length for Q. m = 0.2, and we have "padded" 
each spectrum to 1422 km s _1 by replicating the first 400 km s . Reducing Q, m also 
raises the fluctuation amplitude at z = 3 (since o% is held fixed at z = and there is less 
late time growth for lower Q. m ), but this effect is less important than the effect on the 
Hubble parameter. 

The relative simplicity of the underlying physics and the existence of superb data at 
redshifts that are only sparsely probed by other observables make the Lya forest a po- 
tentially powerful tool for testing cosmological models. Given independent estimates of 
r (e.g., from the quasar luminosity function or the proximity effect), one can use the 
mean opacity of the forest to obtain a lower limit to the cosmic baryon density lfl8L l28ll. 
However, the observational constraints on T are loose, so for higher precision applica- 
tions one must treat T as a free parameter and adjust it to match a single observable, 
usually taken to be the mean opacity. This normalization also absorbs uncertainties in 
Q-b, To, and h. We have followed this practice for each of the simulations in Figure[T]and 
will adopt it in all of our subsequent analyses. After normalizing to the mean opacity, 
the variation about the mean — i.e., the structure of the Lya forest — is a prediction of 
the cosmological model, driven primarily by the structure in the underlying dark matter 
distribution. Ideally, one would like to calculate predictions for each cosmology using 
high resolution hydrodynamic simulations, but for some purposes other numerical or 
analytic approximations may be accurate enough. The adequacy of such approximations 
must be tested on a case-by-case basis, and the improving quality of the observational 
data places ever more stringent demands on the accuracy of the theoretical calculations. 



CONTINUOUS FLUX STATISTICS 

Because of the tight predicted correlation between the local Lya optical depth and 
the local gas density, the most natural statistics for characterizing the Lya forest are 
those that treat each quasar spectrum as a continuous 1 -dimensional field, rather than a 
collection of discrete "lines." We have reviewed this approach elsewhere idQ, and 
here we provide a brief recap with different examples and emphasis. 

Figure |2] shows the power spectrum of continuum normalized flux, converted from 
1-d to 3-d as described by 112211 . and multiplied by k 3 to yield the variance per Ink. The 
curves are noisy because they are based on a single 11.11 1/^ 1 Mpc comoving cube, but 
we use the same phases in each simulation, so the relative behavior should be minimally 
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FIGURE 2. The flux power spectrum. Heavy and light solid curves show the 3-d flux power spectrum 
from, respectively, the SPH simulation and a PM simulation with the same initial conditions. Dotted and 
dashed curves show results with lower and higher matter fluctuation amplitudes, corresponding to (7g = 0.6 
and 1.0 at z = 0. All results are from a single 11.111/z 1 Mpc cube realized with the same Fourier phases. 



affected by noise. At scales k > 0.03 (km s ) , the power spectra turn over because of 
the combined effects of non-linearity, peculiar velocities, and thermal broadening. The 
predictions in this regime are also affected by the finite resolution of the simulations. 
At larger scales, the SPH and N-body+FGPA methods give similar but not identical 
results for the same cosmological model, and the amplitude of the flux power spectrum 
increases with the amplitude of the matter power spectrum, as expected based on the 
FGPA and on Figure [T] 

On large scales, the shape of the 3-d flux power spectrum is similar to that of the 
linear matter power spectrum P m (k) I22U 291. The close connection between the shape 
and amplitude of the flux power spectrum and the shape and amplitude of P m (k) is 
the basis of Croft et al.'s ll22ll method for recovering P m (k) from Lyct forest data. 
Applying this method to a sample of 30 Keck HIRES spectra and 23 Keck LRIS 
spectra yields a matter power spectrum in remarkably (or, perhaps, disappointingly) 
good agreement with the "concordance" ACDM model favored by CMB, supernova, 
weak lensing, and low-z large scale structure data ( HHl ; see ll30ll for an independent 
analysis of the same flux power spectrum and 113 11] for an independent comparison to 
other cosmological constraints). McDonald et al. IB2I1 reach similar conclusions from 
a "forward" comparison of hydrodynamic simulation predictions to the flux power 
spectrum measured from eight HIRES spectra. The Lya forest power spectrum tests 
the ACDM model in a previously unexplored regime of redshift and lengthscale. It 
confirms one of the scenario's key predictions, a linear power spectrum that bends from 
the primeval k n towards k n ~ 4 on small scales. The implied constraints on cosmological 
parameter combinations complement those from the CMB and other data. The Sloan 
Digital Sky Survey (SDSS) has already obtained moderate resolution spectra of several 
thousand high redshift quasars, and these will soon provide measurements of the flux 
power spectrum with much greater precision on large scales. 



0.5 

Pu 
on 

2, 



-0.5 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

D D 

FIGURE 3. PDF of the flux decrement D = 1 e~ T , measured from simulated spectra smoothed with 
a Gaussian of FWHM=8 km s (left) or 100 km s (right). Line types as in Fig. [2] 

Figure |3 shows a different statistic, the probability distribution function (PDF) of the 
flux decrement D = I— F = 1 — e~ z ', for the same set of simulations. The measurements 
in the left hand panel are from simulated spectra smoothed with a Gaussian of FWHM = 
8 km s , comparable to the resolution of Keck HIRES or VLT UVES spectra. The SPH 
and N-body+FGPA predictions agree well for o% = 0.8. Models with higher fluctuation 
amplitude have a broader distribution of densities p and a correspondingly broader 
distribution of flux decrements D, with more saturated and low-opacity pixels and fewer 
pixels of intermediate opacity. The differences in the predicted PDFs for D ~ 0.3 — 0.7 
are AlogP ~ 0. 1 — not enormous, but readily measurable at high statistical significance 
with reasonable observational samples. The fraction of saturated pixels (defined here 
by D > 0.96) increases from 7.2% to 8.3% to 9.1% as cr 8 goes from 0.6 to 0.8 to 1.0. 
There are larger differences in the number of nearly transparent pixels, but continuum 
fitting uncertainties make it difficult to measure the PDF accurately near D = 0. One 
can investigate the scale dependence of matter clustering by smoothing the spectra (or 
by observing them at lower spectral resolution), analogous to studying galaxy counts in 
cells of increasing size. Figure |3] (right) shows results for 100 km s _1 smoothing. The 
predicted PDFs are narrower, since smoothing drives pixel values towards the mean, 
but the dependence on Og is similar. The saturated pixel fraction doubles, from 1.9% to 
3.8%, as as rises from 0.6 to 1. 

For cosmological models with Gaussian initial conditions, the shape of the flux PDF 
depends mainly on the amplitude of mass fluctuations, once one has chosen T to match 
the mean opacity, I33II23I1 . The effective physical scale of this amplitude measurement 
is determined by the spectral smoothing or, for spectra that fully resolve the observed 
absorption features, by a combination of thermal broadening and gas pressure effects. 
Models with non-Gaussian initial conditions predict significantly different flux PDFs 
lE3ll . Observational tests to date show good agreement with models that have Gaussian 
initial conditions and a P m (k) amplitude compatible with ACDM predictions iflil I23I 
l32ll . The uncertainties in these tests are comparable to the model differences in Figure[3] 
Reducing them requires careful attention to the effects of continuum fitting and noise and 
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FIGURE 4. Threshold crossing frequency in spectra with FWHM=8 km s (left) and 100 km s 
(right). Curves show the number of times per unit redshift that the spectrum crosses a threshold of flux 
decrement D as a function of the filling factor, the fraction F(D) of pixels with flux decrement less than 
D. Heavy and light solid curves show the SPH and PM results, respectively, for Q. m — 0.4. Dotted and 
dashed lines show PM results for Q m = 0.3 and 0.2, respectively. 



to the accuracy of the theoretical predictions, but it should be possible to obtain tight 
constraints on the normalization of the matter power spectrum, and some information 
on its shape by studying different smoothing scales. Consistency between results from 
the flux power spectrum and the flux PDF can be a sensitive diagnostic for primordial 
non-Gaussianity. 

FigureH|shows the threshold crossing frequency, the number of times per unit redshift 
that the absorption spectrum crosses a decrement threshold D. Following we plot 
dN c /dz as a function of filling factor, the fraction F(D) of pixels with flux decrement 
less than D, which cleanly separates the information in dN c /dz from the information 
in the PDF and makes the model predictions nearly independent of the photoionization 
rate T. As one would expect from Figure[H the threshold crossing frequency increases as 
Cl m decreases because a lower Hubble parameter H(z) "squeezes" redshift separations 
relative to comoving distances. The threshold crossing frequency drops slightly as the 
amplitude of P,„(k) increases and gravitationally driven merging "smooths" structure 
J3, but the effect is small over the range Og =0.6—1. The crossing frequency also 
drops for redder matter power spectra, since these also lead to smoother structure 
H23I1 . For high spectral resolution, dN c /dz also depends on the gas temperature, which 
determines the level of thermal broadening. 

The great promise of the threshold crossing statistic is its potential for constraining the 
Hubble parameter at high redshift. In spatially flat models with a cosmological constant, 
the ratio H(z) /Hq is determined by Q. m . Alternatively, if Q. m is known independently, the 
ratio H(z) /Hq can constrain the equation of state of dark energy 13 51 13611 . Unfortunately, 
dN c /dz also depends on other cosmological and IGM parameters, and it is numerically 
difficult to predict with high accuracy even when the model is fully specified. The 
difference between the PM and SPH predictions in Figure |4| is comparable to the 
Q. m effects themselves, and the SPH result is still affected by the finite numerical 
resolution. On the observational side, accurate measurement of dN c /dz requires high 



signal-to-noise spectra. Furthermore, while moderate resolution spectra can provide 
useful diagnostics for the shape and amplitude of the matter power spectrum 12311 . the 
H(z) application demands high spectral resolution, since smoothing over a scale that is 
fixed in redshift units erases the sensitivity to H(z) (see Fig.H right). Exploiting dN c /dz 
as a probe of dark energy thus represents a theoretical and observational challenge. 



THE GALAXY PROXIMITY EFFECT 

One can also use the observable correlations between Lyman Break Galaxies (LBGs) 
and the Lya forest to study the environments of high redshift_galaxies. Several groups 
have investigated this issue theoretically using simulations and Adel- 



berger et al. have carried out an observational study using an LBG survey in fields probed 
by seven quasar lines of sight ( 14111 . hereafter ASSP). On large scales — cubes of co- 
moving size ~ 13/^ 1 Mpc — ASSP find a clear correlation of galaxy overdensity with 
Lya flux decrement, the expected signature of galaxy formation in overdense environ- 
ments ff3vL L39ll . However, the observed Lya decrement decreases within A,, ps l/i -1 Mpc 
of LBGs (comoving, redshift-space separation), where the simulations predict that ab- 
sorption should be strongest f: 38ll39il40ll . 

Figure |5] illustrates this conflict. In the upper left panel, triangles show the predicted 
mean decrement in bins of A r , while filled circles show the ASSP data points. The 
predictions come from an SPH simulation of a Tl.TXlhr 1 Mpc cube, and we select 
the 40 galaxies with the highest star formation rates to approximate the magnitude limit 
of ASSP's spectroscopic survey (see [39] for details). Plausible random errors in galaxy 
redshifts can reduce the discrepancy at A r « 1 — 2/i _1 Mpc, but they cannot explain the 
results at the smallest scales, especially the innermost data point at A r = — 0.5h~ l Mpc, 
(D) =0.11 H38LI39LI42I1 . The right column of Figure |5] shows simulated spectra along 
12 lines of sight that pass within 0.5h~ l Mpc of a target galaxy, centered on the galaxy 
redshift space position. All but one of these spectra have D > 0.6 at the galaxy redshift, 
and half have D> 0.9. 

The obvious conclusion is that feedback from the observed galaxies reduces neutral 
hydrogen in their immediate surroundings. Local photoionization by the galaxies' stars 
or active nuclei proves insufficient; this argument can be cast in general terms that seem 
difficult to escape [ 39] . The natural alternative is some form of supernova or AGN-driven 
wind [41, 38^ |40J]. Here (and in 042J]) we have investigated highly idealized "wind" 
models in which we simply eliminate all neutral hydrogen in a sphere of comoving 
radius R w around each target galaxy. Squares, pentagons, and hexagons in Figure[5]show 
results for R w = 0.75, 1.0, and 1.5/z -1 Mpc. Stars show a model where we place winds 
around all 641 resolved galaxies in the simulation volume (instead of the top 40 that 
constitute the "observed" sample) and scale the wind volume in proportion to the galaxy 
baryon mass, with a normalization R w = l/i -1 Mpc around the 40th-ranked galaxy. In 
the right column, dotted curves show spectra for the R w = 1 .5/i 1 Mpc model. 

The most important lesson from Figure |5] is that eliminating neutral gas to a real 
space distance R w does not eliminate absorption within a redshift space separation R w 
(e.g., to AV = ±200 km s _1 in the spectrum plots). Peculiar infall velocities allow gas 
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FIGURE 5. Influence of idealized "winds" on the conditional mean flux decrement and the flux power 
spectrum, based on an SPH simulation of a 22.222/i~' Mpc (comoving) cube. Triangles in the upper left 
panel show the mean flux decrement in pixels that lie at redshift-space separation A r from one of the 40 
brightest galaxies in the cube. Results are averaged over 0,5h Mpc (comoving) bins of A r . Squares, 
pentagons, and hexagons show the effect of removing all neutral hydrogen in spheres of comoving radius 
0.75, 1.0, and 1.5/z Mpc around these galaxies before extracting spectra, and stars show a model with 
sphere volume proportional to galaxy baryon mass. Filled circles show the observational estimates of 
ASSP. Sample spectra on the right illustrate the effect qualitatively. Solid curves show spectra along 12 
lines of sight selected to pass within 0.5h~ 1 Mpc (24") of a target galaxy at redshift space position AV = 0. 
Dotted curves show the corresponding spectra for R w = 1 .5/? 1 Mpc. The lower left panel shows the flux 
power spectrum for the unmodified SPH simulation (solid) and the various "wind" models (as marked). 



at larger distances to produce absorption near the galaxy redshift, making the energetic 
requirements on any wind explanation of the observed Lya deficit much more stringent. 
Only the mass-scaled and R w = 1 .5/? 1 Mpc models come close to matching the ASSP 
data. Winds that fully ionize or perfectly entrain gas to such large distances are not a 
natural outcome of hydrodynamic simulations with stellar feedback II43L 1401. and even 
reaching 1.5/z 1 Mpc in ~ 1 Gyr requires a sustained propagation speed ~ 600 km s _1 . 
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FIGURE 6. Influence of winds on the flux decrement PDF (top, as in Fig.|5Jl and the threshold crossing 
frequency (bottom, as in Fig. [4}. Results are shown for the 22.222/z -1 Mpc SPH simulation with no 
modification (solid) and with the "wind" models illustrated in Fig. [5] 



Given the challenges facing the wind explanation, it is worth considering the alter- 
native possibility that extended Lya emission from the target galaxies is "filling in" 
the corresponding region of the absorption spectrum. Steidel et al. 114411 observed two 
extended Lya "blobs" apparently associated with LBGs, with angular extents ~ 15" 
and AB apparent magnitudes 21.02 and 21.14 in the Lya band. Cooling radiation from 
gas settling into massive galaxies at z = 3 naturally produces Lya fluxes of this order 
. The three quasars that contribute to ASSP's innermost data point have G-band 
AB magnitudes of 20.1, 21.6, and 23.4, so if all of a target galaxy's Lya cooling ra- 
diation went down the slit it could potentially replace the quasar flux absorbed by the 
surrounding IGM. However, a 1.4" slit at an angular separation AO ~ 15 — 20" from a 
galaxy should intercept at most ~ 1.4/2A0 ~ 0.03 — 0.05 of the galaxy's extended Lya 
flux, at least on average, so this explanation seems to fail by 1 — 2 orders of magnitude. 
Furthermore, a fourth pair involving the G = 17.8 quasar Q0302-0019, shows no sign of 
absorption near the galaxy redshift, and in this case the quasar is clearly too bright for 
galaxy emission to compete with it. (This pair and two others are dropped from ASSP's 
(D) calculation because of possible Ly/3 contamination.) At this point, the Lya emis- 



sion explanation seems unlikely; it can be conclusively ruled out by observing more 
close pairs involving bright quasars or by obtaining symmetrically placed spectra away 
from the observed quasars to search for galaxy emission. 

Assuming for now that winds are the correct explanation for the observed Lya deficit, 
one can ask whether they completely spoil the picture painted in §§1 and 2. The lower 
left panel of Figure |5] shows the flux power spectrum for the various "wind" models. 
Eliminating neutral hydrogen to distances R w = 0.75 or l/i -1 Mpc around bright galaxies 
has only a small impact on the flux power spectrum, because the filling factor of the 
"bubbles" is small and absorption close to the galaxies remains nearly saturated in any 
case. However, these models also do not explain the ASSP results. In the two more 
extreme models, winds suppress the flux power spectrum on large scales by 0.1 — 0.2 
dex, comparable to the 0.1 -dex la uncertainty that [29] quote for the normalization 
of the matter power spectrum. Thus, winds of this magnitude could have systematic 
effects on P m (k) determinations that are significant relative to the present observational 
uncertainties. The influence of winds on the flux PDF and threshold crossing frequency 
is smaller, as shown in Figure|6l Even the more extreme models have a negligible impact 
on dNc/dz, and they only slightly alter the shape of the PDF. The most significant effect 
is on the fraction of saturated pixels for 100 km s~ ! smoothing, which drops from 2.5% 
for the no-wind case to 1 .9% for R w = 1 .5/i _1 Mpc and 1 .3% for the mass-scaled model. 



PROSPECTS AND CHALLENGES 

From the above discussion, it is clear that one immediate challenge is to better under- 
stand the galaxy proximity effect. Figures |5] and ]§\ show that winds can reach substantial 
distances from bright galaxies without having much impact on the global statistics of 
the Lyoc forest (see also [4-3, 40]), but even the more extreme wind models considered 
in §3 do not reproduce the ASSP results very well. Simulations with more realistic wind 
physics IB8l |43l 1401 can shed further light on this problem, but major progress will have 
to await further observational studies, since current inferences rest crucially on a handful 
of galaxy-quasar pairs. 

Studies with Keck HIRES and VLT UVES have provided a (still growing) trove 
of high quality data on the Lya forest at z ~ 2 — 4. The fluctuating IGM perspective, 
furthermore, shows that large samples of moderate resolution spectra can be a powerful 
resource for studying structure on large scales, since continuous flux statistics do not 
require resolution of individual "lines." Such samples can be assembled quickly with 
large telescopes, and the SDSS is producing an enormous sample at resolution R ~ 2000 
in the course of its normal operations. One of the observational frontiers is the use of 
correlations across multiple lines of sight. Quasar pair studies provided the first decisive 
evidence for large coherence scales of absorbing structures JEIHIUBI], but larger samples 
allow more ambitiou s go als, such as using the Alcock-Pacyznski test |Et1i to measure 
spacetime geometry ||48|j49|,[50], improving measurements of the flux power spectrum 
with cross correlations 15 ill , and mapping large scale 3-dimensional structure at high 
redshift JHB 

Analyses of existing data have already led to an important cosmological conclusion, 



namely that models with matter clustering similar to that of "concordance" ACDM at 
Z ~ 3 are consistent with the observed Lya forest while models with substantially differ- 



ent clustering amplitudes or suppression of small scale power are not [18, 541 I55LI23L I32, 



SEilSEjJIlIl- They have also provided constraints on the temperature of the dif- 
fuse IGM [59[ 6C JfH] and indications of helium reionization at z ~ 3.2 ll6lll62Ll63l l64ll. 
However, to do justice to the quality and quantity of data and keep pace with the tighten- 
ing observational constraints from other observables, we must play for higher stakes. It 
looks possible in principle to achieve precision of 5 — 10% on quantities like the matter 
fluctuation amplitude and H(z), but even without the potential complications of galaxy 
feedback, inferences at this level require more extensive theoretical modeling. Many ef- 
fects that are unimportant at the 25% level — differences between approximate methods 
and full hydrodynamics, numerical resolution and box size limitations, spatial fluctua- 
tions in IGM temperature, details of continuum determination — may become critical 
at the 5 — 10% level. Despite these challenges, the Lya forest is the most promising 
tool we have for precision cosmological measurements at z ~ 2 — 4. These measure- 
ments might, in the end, simply confirm the cosmological model favored by other data, 
but complementary constraints have the potential to break parameter degeneracies and 
thereby reveal subtle quantitative discrepancies. These in turn could yield insight into the 
nature of dark energy, the mechanisms of inflation, or some other fundamental aspect of 
our universe. 
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